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Introduction 

l,3.5-trinitro-l,3,5'triazacyclohexane  (1;  RDX)  and  l,3,5,7-tetranitro-l,3.5,7-tetraaza- 
cyclooctane  (2;  HMX)  are  highly  energetic  cyclic  nitramines.*’^  Both  can  be  represented  by  the 
formula  (CH2NNC)2)n-  The  smallest  member  of  this  series,  l,3-dinitro-l,3-diazacyclobutane  (3), 
has,  to  our  knowledge,  not  yet  been  synthesized.  A  possibly  significant  feature  of  3  that  is  not 
present  in  1  or  2  is  strain  due  to  the  four-membered  ring.  This  should  increase  the  molecule's 
energy  content,  although  the  effect  is  likely  to  be  reduced  by  a-conjugation  of  the  aza  nitrogens' 
lone  pairs.^ 


0,N— N 
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H,C— N 

?■  \ 


N— CH, 


N—NO, 


N - CH, 
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A  useful  measure  of  energy  content  is  the  heat  of  formation,  and  one  of  our  present 
objectives  was  to  compute  this  for  3.  We  also  extend  earlier  preliminary  studies  of  two  possible 
decomposition  steps:^’^  (a)  homolytic  rupture  of  an  N-NOa  bond,  and  (b)  dissociation  of  the  ring 


into  two  H2C=N-N02  fragments,  passing  through  the  transition  state  4. 

Methods 

HF/6-31G*  optimized  geometries  and  vibration  frequencies  were  computed  for  3  -  7  with 
the  GAUSSIAN  92  program;®  the  frequencies  were  scaled  by  1/1.12.®’^  (For  5  and  6,  the 
unrestricted  Hartree-Fock  approach  was  used.)  Single  point  runs  were  then  carried  out  for  each 
structure  with  the  non-local  density  functional  (DF)  code  deMon,®  using  the  Gaussian  DZVPP 
basis  set;  the  exchange-correlation  energies  were  calculated  by  means  of  the  generalized  gradient 
approximation  (GGA).^’*®  These  DF  results  were  used  in  computing  the  reaction  energetics. 

Standard  gaseous  heats  of  formation  were  obtained  by  computing  AE  (DF)  for  the 
formation  of  each  molecule  from  its  elements,  combining  this  with  translational,  rotational  and 
vibrational  contributions  and  converting  the  result  to  at  25°C  using  the  ideal  gas 
approximation.* Empirical  bond  and  group  correction  terms  were  also  included.* 

Results 

The  HF/6-31G*  optimized  geometries  of  3  -  7  are  given  in  Table  I.  Although  experimental 
values  are  known  only  for  NO2  (6),  data  relevant  to  3  and  5  are  in  crystallographic*^’*^  and 
electron  diffraction*^  studies  of  RDX  (1)  and  1,3,3-trinitroazetidine,  and  an  MCSCF  structure  is 
available  for  7.*®  In  general,  there  is  good  agreement,  although  the  HF/6-31G*  procedure  tends  to 
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underestimate  N-0  bond  lengths  by  roughly  0.03  A.  Our  structures  for  3  and  5  -  7  do  correspond 
to  local  energy  minima  (no  imaginary  firequencies)  and  4  is  a  transition  state  (one  imaginary 
frequency).  The  HF/6-31G*  and  the  DF-GGA/DZVPP  energies  of  3  -  7  are  in  Table  H. 

The  strain  energy  of  3  was  estimated  by  the  isodesmic  reaction  procedure,^  which  involves 
writing  a  hypothetical  chemical  process  in  which  the  number  of  bonds  of  each  type  remains  the 
same  but  their  mutual  relationships  have  changed.  AE  for  such  a  reaction  is  interpreted  as 
revealing  any  deviations  from  bond  energy  additivity,  due  to  strain,  resonance  stabilization,  etc. 

At  the  HF/6-31G*  level,  taking  H2N-NO2  and  H3C-NH2  to  be  the  reference  molecules  and 
using  energies  obtained  in  this  work  and  from  archives,^^  we  find  the  strain  energy  of  3  to  be  8.6 
kcal/mole.  This  is  in  good  agreement  with  our  earlier  HF/3-21G  value  of  9.5  kcal/mole.^  As  a  test, 
we  also  computed  the  HF/6-31G*  strain  energy  for  cyclopropane,  using  an  analogous  isodesmic 
reaction,  and  obtained  26.3  kcal/mole,  which  is  close  to  the  quoted  value  of  28.3  kcal/mole. 

Our  estimated  standard  gaseous  heat  of  formation  at  25°  C  for  3  is  -1-52.0  kcal/mole.  We 

believe  this  to  be  reliable,  since  the  very  similar  RDX  (1)  was  part  of  our  data  base  in  determining 
the  correction  terms,  and  our  procedure  reproduces  its  measured  AHf  to  within  0. 1 
kcal/mole.  Our  AHf  for  3  exceeds  that  of  RDX  by  6.2  kcal/mole;  this  may  reflect  the  strain  in  3. 

Using  the  data  in  Table  H,  with  zero-point  energies  taken  into  account,  the  DF  N-NO2 
dissociation  energy  is  36.6  kcal/mole.  The  overall  AE  for  the  ring  fragmentation  is  17.0  kcal/mole, 
but  the  activation  barrier  for  this  process  is  44.2  kcal/mole.  N-NO2  dissociation  and  ring 
fragmentation  are  therefore  roughly  similar  in  terms  of  their  energetic  requirements. 

Discussion 

The  heat  of  formation,  expressed  on  a  mass  basis,  is  a  key  determinant  of  the  detonation 
and  propellant  properties  of  an  energetic  material.^  Thus  it  is  noteworthy  that  our  value  for  3, 

351  cal/g,  is  70%  larger  than  that  of  RDX  (1),  206  cal/g. 

The  N-NO2  dissociation  energy,  36.6  kcal/mole,  is  somewhat  less  than  the  40-50  kcal/mole 
range  typical  of  N-NO2  bonds,^*'^  possibly  reflecting  some  release  of  strain  upon  breaking  the 
bond.  In  view  of  evidence  that  rupture  of  the  N-NO2  bond  is  an  important  factor  in  nitramine 
decomposition,^* the  weakness  of  these  bonds  in  3  suggests  that  it  might  possess  a  rather  high 
sensitivity  to  shock  and/or  impact.  A  similar  conclusion  is  reached  on  the  basis  of  a  correlation 
relating  nitramine  shock  sensitivity  to  the  N-NO2  bond  lengths  and  the  molecular  weight;  ^  3  is 
predicted  to  have  a  shock  sensitivity  very  close  to  those  of  RDX  (1)  and  HMX  (2).  Nevertheless, 
its  relatively  high  energy  content,  as  reflected  in  its  heat  of  formation,  should  stimulate  continued 
synthetic  efforts. 
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Table  I.  HF/6-31G*  optimized  geometries. 


Molecule 

Distances  (A) 

Angles  (deg) 

OjN 

C-N 

1.466 

N-C-N 

87.3 

^  N 

N - CH 

N-N 

1.355 

C-N-C 

92.7 

1  1 

N-O 

1.191 

C-N-N 

119.9 

1  1 

3 

N-N-O 

116.4 

HjC 

O-N-0 

127.2 

NO2 

N-C-N-N 

127.1 

C-N-N-0 

34.6,  147.7 

C-N 

1.343 

N-C-N: 

92.0 

N 

1 ---  •  ^ 

C-N 

1.991 

C-N-C: 

88.0 

1  1 

1 _ 1 

4 

N-N 

1.369 

C-N-N:' 

114.1 

HjU 

N-Oa 

1.191 

N-N-Oa.*’  115.0 

NO2 

N-Ob 

1.197 

N-N-Ob:*’  119-4 

O-N-O: 

125.6 

N-C-N-C: 

0.0 

C-N-N-Oa:‘»  156.7 

C-N-N-Ob:'’  24.2 

C-N(N02) 

1.466 

N-C-N 

90.3 

, 

C-N 

1.459 

C-N(N02)-C 

88.8 

N 

1  1  ^ 

N-N 

1.349 

C-N-C 

89.3 

1  1 

N-O 

1.193 

C-N-N 

121.2 

HjC— N 

N-N-O 

116.6 

NO2 

O-N-O 

126.8 

C-N-C-N 

8.5 

N-N-C-N 

134.9 

C-N-N-O 

35.9,  145.8 

•NO2 

6 

N-0 

1.165 

O-N-O 

136.1 

(1.197)a 

(133.8)a 

C-N 

1.254 

C-N-N: 

113.7 

N-N 

1.417 

O-N-O: 

126.6 

H2C=N-N02 

7 

N-Oa 

1.180 

N-N-Oa: 

113.9 

N-Ob 

1.193 

N-N-Ob: 

119.5 

C-N-N-Oa: 

157.1 

C-N-N-Ob: 

24.1 

^Experimental  values  for  NO2.  in  parentheses,  from  Harmony,  M.  D.;  Laurie,  V.  W.; 
Kuczkowski,  R.  L.;  Schwendeman,  R.  H.;  R^amsay,  D.  A.;  I^vas,  F.  J.;  Lafferty,  W.  J.;  Maki, 
A.  G.  J.  Phys.  Chem.  Ref.  Data,  1979, 8,  619. 

'’These  angles  are  within  each  H2C=N-N02  portion. 


Table  n.  Calculated  energies. 


Molecule 

Total  Energy* 

Zero-point  energy 

(HF/6-31G*), 

(kcal/mole) 

HF/6-31G* 

DF-GGA/DZVPP//HF/6-3 IG* 

(haittees) 

3 

-594.97670 

-599.02944 

58.13 

4 

-594.86945 

-598.95329 

54.53 

5 

-390.91306 

-393.62772 

47.53 

6 

-204.0.3149 

-205.33525 

5.52 

7 

-297.47560 

-299.49753 

26.76 

^This  corresponds  to  energy  minimum  at  0°K. 


